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ABSTRACT 


, i j is given of performance characteristics of higher oraer 
; of Runge-Kutta type, performance predictors for the tznc 
;e machine and for the size of the error are developed. These 
:• not supposed to give precise information, but supporting 
•_ in the note that indicates useful information is obtained 
Lctors. The predictors and the data indicate that the usual 
should be Shanks 1 , formulas of order six, seven, and eight. 
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IIvTRODUCTICX. 


The -writer of this note has developed higher order formulas for the 
r .".proximate solution of ordinary differential equations. The previous 
formulas of this type extended only through the sixth order (tvo of the 
sixth order vere given hy Huta, each .requiring eight evaluations of the 
function;). Among the new ones developed by the -writer, programmed by 
Mr. Pred Calhoun, tested, and available for use in the Marshall Space 
Flight Center Computation Division are several of the sixth order recuir 

:nd one of the eight order 
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There are three important factors involved an the use of one of the 
formulas; namely, the formula itself, the differential equation, and the 
computing machine. In this classification, the programming has been in- 
cluded in the first factor. The pertinent characteristics of the fermul 
are the time r (in minutes throughout this note) of one- srep, rhe order 
of the formula, and the number e of evaluations of the. function used in 
defining the differential equation (for convenience, it is ass: .od that 
the differential equation is expressed in the form y’ = f(x,y)_ / . The 
characteristics of the machine are the rime t a of one "ideal" arirh_-.e-.ic 
op erariom ami rhe rime t e cf one evaiuacion of rhe function. (r g cbviousl 
depends or. rhe differential equation also).. 


Let h denote.- the step size, s the number of seeps involved in t, ; .u 
.. elution procedure, T « ts. the total time, and E » es the total run' or 
evaluations. 

For the time of one step, one deduces immediately the (approximate) 
equation 

t «« et e + e(e+l)t a , 

since e(e~4) is the (approximate) number of arithmetic operations require: 
in the formula (this does not take into account the presence of zero co- 
efficients in the formula). For the total time, one deduces the (approxi- 
mate) equation 


T = 2t e + E(e+4)1 


( 1 ) 


lie quantity t a is (approximately) constant for s 
liven type of precision programmed on the machin: 
vary significantly for different differential o _ 
Equation (1) leads to the following equiv_l: 


given machine and a 
ihe quantity .t 3 will 
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T2 = E 2 


T-, / I 

+ ( e a- e i. -aj > 


( 2 ) 


where It, e a refer to one formula, T 2 , S 2 , e 2 refer to another for- 
mula and the same differential equation is involved. It is assumed in 
vhas follows that the differential equation is fixed, unless the contrary 
is explicitly stated. 

It is seen from equation (2) that t a may he calculated empirically 
from the results of two runs on the machine of two formulas' with distinct 
e numbers. Ihis has teen done on the IBM 7090 in extended"precision and 
t a is (approximately) (kp)10“ . For a particular fixed differential equa- 
tion, one run on the machine with a part_eular formula detema.-es IS, 2^,0 
Then equation (2) can he used to predict the .time T required 00 m.ake a run 
with E evaluations for other formulas. . ’ 
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PREDICTED ACTUAL 
• TBS TBS 

PREDICTED 

ACTUAL 

FORMULA. 

TIPS (mbsutes) (miiutes) 

ERROR 

ERROR h E 

1. RUNGE-XUTTA 

(M) 

11.66 9.81 

1C " 7 

.10 io" s .oi 72000 

.2. RUIvC-2-KUTTA 

(M) 

5.89 ’ 3.27 

io “ 3 

• pi 10 • op 

3. RUT J L-S - ICUT’TA 

(M) 

1.94 1.65 

10" 4 

.51 10~ 3 .C6 12GC0 

SHARKS 

(o,7)#l 

l . 1 7 l • z>z> 

io ‘ 4 

.17 10“'*'. 15 Sloe 

^ c* T~rc 

( 6 , 7)#2 

1.1.7 i.Lg 

- 0 -* 
— .w 

.16 10" 3 .13 84CC 

a. £ EASES 

(6,7^3 

1. . ' 


.cv ii"-.-_= 

?• SHARKS 

(°,7)#4 

J m J • GO 

— V 

-■ “> r\“" 79 ~ " ' *'■ 

• — ^ __ W » s/v/ w — 

8 . SHARKS 

(7,9)#1 

2.49 . 2.53 

1 c ' 3 

.15 io" s .i 2 13500 

9- SSsSBS 

(7,9)#2 

2.^9 2.52 

io- s 

.31 1-“ 5 .12 15500 

10. SHARKS 

(8,12) 

4-75 >-73 

io“ 7 

.11 1 ’'.09 24000 • 

The predi 

cted error is taken to he 10' 

-n . ■ . 

, where h 

3 =(.d 1 d 2 ...)lO“ r * 

and 0 < di< 9* 
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the various 
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following table. 
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(M) 


(5,3) 


.14 10“~ 

(o,7) 


(M) 


12 10“ 3 

(.7,9) 


(6,7) 


.20 10“ = 

(8,12) 


(7,9) 


.20 1C“ S 

(5,6) 


(L,L) 

# 

- l *4- 

f £ rr\ 
V°, ! 1 


(5,6) 


.29 1 C~“ 

(7,9) 


(^) 


. 48 10“ 


» 02T ■ ■— »- ^ r] *1 , Z* wrc' M WMk.« cCmUu.uH 


• W»J* Z,Z*L> -u. — -- 


y" , y - A., x c 


■ ^y 2 . xy 

equation (2) becomes (approximately) 


= y 0 » y'o = 1 


• T = — (e + 32) . 
10 7 


( 5 ) 


Suppose new that the sane interval as well as the sane differential 
-cuation is used. Without loss of generality, it nay he assumed than 
S > Qt. and e > e^ and that the (0,e) formula is run for only one step. 

0J--1 h 0 -*-~ 

Suppose further that h‘ 2 = (j) 1,2 , that is, that the errors in the 
two formulas are (approximately) equal. Then 

Qi-6 

s-, = h - 2 . 
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if and only if 
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■0 1 


h 8l+ 2 


= Si > ~ 
- Z 1 


is inequality, one easily deduces 
he used when the error desired is to he , 
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.nan the critical value 
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e 1 + (e + 4)t r 1 


(M 


where t r is defined hy the equation t r t e = t a . 

The equation (2) and the 'critical value (h) are the-usefu_ informa- 
tion referred to in the introduction. Ke. cher should he interpreted as 
giving precise information. However, as s. own hy the data' given below, in 
cases where the differential equation involves a function that reason 
.sly well- behaved, 11 the predictions may he expected to he reasonably ac- 

4. sate, moreover, it is reasonable to expect that neither of the: o for 

.veins will be good approximations when • 


nr edict ions are nor. 


tne 
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?rcm onese predictions, it follows that there is no need f.feh 
order formulas since the writer has proved that a fifth order formula 
r.-.rct require at least six evaluations. It can also he seen that chi- pr_- 
» eren ce lisu is biased toward the lover order formulas hy examining rur.s 
1 ar.d 7; runs 3 ur.d 6; runs 7 and 9 . 

■These results clearly indicate that (l) Hystrcm 1 s formula and Plata's 
formulas are out-dated and should not he used, (2) the Hunge-Kutta formula 
snculd he used only when (relatively) rough approximations are desired, 
vt) tne usual formulas used should he Shanks 1 formulas of order six, 
seven, and eight. 

These results are shown graphically hy the following diagram. 
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. 12 ' 10 “° . 1 , 10 *- 

(M) .. (5,5) 
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accuracy 
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4. Sill' IF: Y 


1. Definition of Symbols. 

h step size used in the solution 

s the number of steps 

t a • machine time for one "ideal" arithmetical operation 

machine time for one "ideal" evaluation of the function 
the number of evaluations in a particular formula 
the order of approximation of a particular formula 
a formula of order Q with e evaluations 
the number of minutes required on machine ' for one step 
T = ts is the total time for s steps (time in minutes ) 
2 = es is the total number of aval 
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2. Fundamental relations (approximate) . 
r* _ 

Ta n-Sa I rr* *r (®2 - e^to, j 
13 J 

t a «• (^5)10 for extended precision on IBM 7090 
I = 10~ a , .where h Q ~ (‘d x d z . . . )10 _ri and 0 < d x < 9* 


Critical Value ' = 



L 1 4- 1 

(ex + K 

it r ] 

[i + 

(e + ^)t r ] ( 


( 29 j. 1)(28 1 -r 1 ) 
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Use formula (Sj e) when error desired is less than critical value. 
Critical .value for (6,7) and (h_,h) formulas is approximately 
1.2 10 “ 4 , 
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